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ABSTRACT 

We present a detailed study of the metal-enriched circum-galactic medium of a massive galaxy at 
z = 3 using results from the "Eris" suite of new cosmological hydrodynamic "zoom-in" simulations 
in which a close analog of a Milky Way disk galaxy arises at the present epoch. The reference run 
adopts a blastwave scheme for supernova feedback that generates galactic outflows without explicit 
wind particles, a star formation recipe based on a high gas density threshold, and metal-dependent 
radiative cooling. Eris main progenitor at 2; = 3 resembles a "Lyman break" galaxy of total mass 
Mvir = 2.4 X 10^^ M0, virial radius i?vir = 48 kpc, and star formation rate 18 M© yr~^, and its metal- 
enriched circum-galactic medium extends as far as 200 (physical) kpc from its center. Approximately 
41%, 9%, and 50% of all gas-phase metals at z = 3 are locked in a hot (T > 3 x 10^ K), warm 
(3 X lO^K > T > 3 X 10^ K), and cold (T < 3 x 10"^ K) medium, respectively. We identify three 
sources of heavy elements: 1) the main host, responsible for 60% of all the metals found within 3i?vir ; 
2) its satellite progenitors, which shed their metals before and during infall, and are responsible for 
28% of all the metals within 3i?vir , and for only 5% of those beyond 3i?vir ; and its satellite dwarf 
companions, which give origin to 12% of all the metals within 3i?vir and 95% of those beyond 3i?vir • 
Late (2; < 5) galactic "superwinds" - the result of recent star formation in Eris - account for only 
9% of all the metals observed beyond 2i?vir , the bulk having been released at redshifts 5 < 2; < 8 by 
early star formation and outflows. In the IGM, lower overdensities are typically enriched by 'older', 
colder metals. Heavy elements are accreted onto Eris along filaments via low-metallicity cold inflows, 
and are ejected hot via galactic outflows at a few hundred kms~^. The outflow mass- loading factor 
as a function of redshift ranges between 0.1 and 1.2, and shows no correlation with the mass of the 
host over the interval 5.0 x 10^ M© < Mvir ^ 2.4 x 10^^ M©. 

Subject headings: galaxies: evolution - galaxies: high-redshift - intergalactic medium - method: nu- 
merical 



1. INTRODUCTION 

Studies of the ionization, thermodynamic, and kine- 
matic state of heavy elements in circum-galactic gas hold 
clues to understanding the exchange of mass, metals, 
and energy between galaxies and their surroundings. 
The distribution of observed metals in the intergalactic 
medium (IGM) is highly inhomogeneous, with a global 
cosmic abundance at z = 3 of [C/H] = —2. 8 ±0.13 for gas 
with overdensities 0.3 < S < 100 (ISchave et an [2003). 
The characteristic epoch of this enrichment and its main 
donors remain uncertain. Late supe rnova-driven "super- 
' winds" from massive g alaxies (e.g., I Aguirre et al.|[2QQTI : 
I lAdelb erger et al."2003K early outflows from dwarf galax- 
iesje.g., Dekel & Silk 1986; MacLow & Ferrara 1999; 
Madau. Ferrara. Reesl 120011: IMori. Ferrara. &: Mada 



and intergalactic gas, but their relative contributions are 
not constrained by present data, and involve complex 
baryonic processes that are difficult to model. 

Observations clearly show that galactic-scale outflows 
with velocities of several hundred kms~^ are ubiq- 
uitous in massive star-forming galaxies at high red- 
shift and in starburst galaxies in the local universe 
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Furlanetto fc Loeb 2003), quasar-driven winds (e.g. 



Scannapieco & Oh 2004), and the ejection of gas during 
the merging of protogalaxies (Gnedin 1998) are all likely 
to have left some chemical imprint on circum-galactic 
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(e.g^lHeckman, Lee, & Milev'"1990'; 'Pettin i et all 
Martin 2005; Veilleux, Cecil & Bland-Haw thornI 
Weiner et al. 2009; Steidel et al.i i2010» ). They also in- 
dicate that the metal content of the z 3 IGM is closely 
correlated with the positions of nearby galaxies: O vi 
systems are found to be strongly associated with known 
Lyman-break galaxies (LBGs), most detectable C iv sys- 
tems {New ^ 10^^"^ cm~^) lie within 1 proper Mpc 
from an LBG, and roughly one- third of all "intergalactic" 
absorption lines with New ^ 10^^ cm~^ are produced 
by gas that lies within ^80 p roper kpc from an LBGs 
(jAdelberger et al.l 120031 l2005f i. Somewhat puzzling, the 
comoving mass density of C iv ions i n the IGM is found 
to drop by a facto r of 4 at z > 5 (|Rvan- Weber et al] 
^2009"; 'Simcoe etaU[20n). 

Understanding the galaxy-IGM ecosystem, i.e. the 
role of inflows, outflows, star formation and active galac- 
tic nucleus (AGN) "feedback" in governing the gaseous 
and metal content of galaxies and their environment, 
is the goal of many recent theoretical efforts. Hydro- 
dynamical simulations of galaxy formation over large 
cosmological volumes, while advancing rapidly over the 
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past few years (e .^. Cen, Nagamine , fc Ostriker 
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OpDenheimer fc D ave 2008, 2009; Wi ersma et all 120091 : 
Shen et all 120101 : ICen fc Chisari .20111 typically suffer 
from poor spatial and mass resolution - which limits their 
ability to follow self-consistently the venting of metals by 
small galaxies and the transport of heavy elements from 
their production sites into the environment. And while 
with suitable wind prescriptions these simulations have 
shown that galactic outflows can potentially enrich the 
IGM to approximately the observed levels, it is not clear 
whether the same prescriptions are also able to produce 
realistically looking galaxies. To shed some light on the 
objects and processes responsible for seeding the circum- 
galactic and intergalactic medium with nuclear waste, 
we follow here a different approach. We present results 
from a new cosmo logical A^-body /smooth particle hydro- 
dynamic (SPH) suite of simulations of extreme dynamic 
range in which, for the first time, a close analog of a 
Milky Way disk gal axy has been shown to arise at z = 
(jGuedes et al.ll2011f ). Termed "Eris", the simulations fol- 
low the assembly of a massive galaxy halo with a total of 
nearly 20 million particles (gas + dark matter + stars) 
within the final virial radius, and a force resolution of 
120 pc. The feedback from an active galactic nucleus is 
neglected, and a star formation recipe is adopted based 
on a high gas density threshold. Heating by supernovae 
occurs in a clustered fashion, and the resulting pressure- 
driven outflows at high redshifts remove low angular mo- 
mentum metal-enriched gas. As shown below, our de- 
tailed study of the circum-galactic medium of Eris' main 
progenitor at z = 3 reveals that late galactic superwinds 
- the result of recent star formation - account for only a 
small fraction of all the metals found between ~ 100 and 
200 (proper) kpc from the center of the main host, and 
resolves two other sources for the heavy elements found 
in Eris' environment: its satellite progenitors - which 
deposit their metals before and during infall - and its 
orbiting dwarf companions. 

The paper is organized as follows. In § 2 we describe 
the Eris simulations. § 3 presents a detailed study of Eris' 
circum-galactic medium. The origin of circum-galactic 
metals and the role played by satellite progenitors and 
companions in polluting Eris' environment are discussed 
in § 4. The properties of the supernova-driven outflows 
from the main host are analyzed in § 5. Finally, in § 6 
we summarize our main results. Throughout this work, 
we use the metal mass fraction Z c^) = 0.0142 for the solar 
abundance (jAsplund et al.ll2009f ). 

2. THE ERIS SIMULATIONS 

The Eris suite of simulations was performed in a 
Wilkinson Microwave Anisotropy Probe 3-year cosmol- 
ogy running the p arallel TreeSPH code Gasoline 
(IWa dslev et al.l I2004D. Details of the main simulation 
are given in fGuedes et al.l ()2011l ) , and are briefly summa- 
rized here. The high-resolution region, 1 Mpc on a side 
at z = 0, is embedded in a low-resolution, dark matter- 
only, periodic box of 90 comoving Mpc on a side, and 
contains 13 million dark matter particles and an equal 
number of gas particles, for a final dark and gas particle 
mass of moM = 9.8 x 10^ Mq and mspH = 2 x 10^ M©, 
respectively. The gravitational softening length, e^, was 
fixed to 120 physical pc for all particle species from z = 9 
to the present, and evolved as 1/(1 -\- z) from z = 9 



Fig. 1. — The optical/UV stellar properties of Eris at z = 3. 
The images w ere created with the radiative transfer code Sunrise 
f Jonsson' '2006), and show a rest-frame and NUV stellar 

composite of the simulated galaxy seen face-on {top panel) and 
and edge-on (bottom panel). A Kroupa IMF was assumed. The 10 
kpc bar in the bottom panel is in proper units. 

to the starting redshift of z = 90. Compton cooling, 
atomic cooling, and metallicity-dependent radiative cool- 
ing at low temperatures (Mashchenko et al. 26oEt) are 
included. A uniform UV background modifies the ioniza- 
tion and excitation state of th e gas and is implemented 
using a modified version of the iHaardt Madaul (jl996r ) 
spectrum. 

Star formation occurs stocastically when cold (T < 
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3 X 10^ K), virialized gas reaches a threshold density 
nsF = 5 atoms cm~^ and is part of a converging flow. It 
proceeds at a rate 

dp^/dt = O.lpgas/^dyn CX Pg^g (1) 

(i.e. locally enforcing a Schmidt law), where and pgas 
are the stellar and gas densities, and tdyn is the local 
dynamical time. Each star particle has an initial mass 
m* = 6 X 10^ Mq and represents a simple stellar popula- 
tion with its own age, metallicity, and a iKroupa et al.l 
(1993) initial stellar mass function (IMF). Star parti- 
cles inject energy, mass, and metals back into the ISM 
through Type la and Type II SNe and ste llar winds, fol- 
lowing the recipes mson eTall (|2QQ6l ). Each SN II 
deposits metals and a net energy of 0.8 x 10^^ ergs into 
the nearest neighbor gas particles, and the heated gas 
has its cooling shut off (to model the effect of feedback 
at unresolved scales) until the end of the momentum- 
conserving phase of the SN blastwave, which is set by 
the local gas density and ten iperature and by the to - 
tal amount of energy injected (iMcKee fc Ostrikerlll977[ ). 
Cooling is turned off only for those particles within the 
blast radius (for a maximum of 32 neighboring gas par- 
ticles), and no kinetic energy is explicitly assigned to 
them. For the typical conditions of star-forming clouds 
resolved in this study, this translates into regions of size 
Re ^ 30£^5i^^ pc heated by individual SNe and having 
their cooling shut off for a timescale t^; ~ 5 x l{)^Er^i^ yr, 
where E^i = £^/10^^ergs. The energy injected by many 
SNe adds up to create larger hot bubbles and longer shut- 
off times. The advantage of this physically- motivated 
feedb ack model compared to oth er ^^sub-^rid" schemes 
(e.g. ISpringel fc HernquistI 120031 : lOppenheimer fc Pavel 
l2008f ) is that it keeps galactic outflows hydrodynam- 
ically coupled to the energy injection by SNe (albeit 
with a delay). In combination with a high gas density 
threshold for star formation (which enables energy de- 
position by SNe within small volumes), this scheme has 
been found to be key in pr oducing realistic dwarf galax- 
ies ( Governato et al. |[2010f ) and late- type massive spirals 
(jGued es et al . 2QlII). 

Metal e nrichment from SN II and SN la follows the 
model of iRaiteri et al.l (| 19961 ). For SN II, metals are 
released as the main sequence progenitors die and dis- 
tributed to gas within the blastwave radius. Iron and 
oxy gen are produced a ccording to the following fits to 
the lWooslev fc Weaverl (|1995D yields: 

/ X 1.864 

Mpe = 2.802 X 10-M^J Mo, (2) 

and 

/ \ 2.721 

Mo = 4.586 X 10-"^ f ^ j ^0- 

Each SN la produces 0.63 M^;^ of iron and 0.13 Mq of 
oxygen (Thielemann et al. 1986) and the metals are dis- 
tributed between the nearest gas particles. Radiative 
cooling is not disabled following a SN la. Stellar wind 
feedback was implemented based on iKennicutt et al.l 
(|1994f ). a nd the retur ned mass fraction was determined 
following IWeidemannI ([l987 ). The returned gas has the 
same metallicity as the star particle. 




Fig. 2. — Projected gas density (top panel) and gas particle metal- 
licity (bottom panel) of Eris's circum-galactic medium at z = 3 in 
a cube of 500 (proper) kpc on a side. The galaxy's stellar disk is 
edge-on in this projection. 

The reference simulation discussed here is a new 
twin Eris run ("ErisMC") that includes metallicity- 
depen dent radiati v e coo ling at high temperatures fol- 
lowing (She^^eT^l] ()2010f ). Aside from high-temperature 
metal cooling, Eris and ErisMC were run with exactly 
the same setup. For convenience, we will refer to ErisMC 
simply as Eris in the following. 

3. ERIS' CIRCUM-GALACTIC MEDIUM 

Figure [1] shows the optical/UV stellar properties of 
Eris at z = 3, just after its last major merger. The 
mock images were created using the radiation transfer 
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Fig. 3. — Same as Fig. [21 showing the mass- weighted gas metaUicity and velocity field in a 500 x 500 x 10 kpc^ slice. The galaxy center 
is indicated by the plus sign at coordinates = (0,0), and its virial radius is marked by the black circle. The color bar indicates the 

mean gas metallicity within the slice, with pristine gas shown in gray. The arrows depict the mass-weighted peculiar velocity vectors of the 
gas relative to the main host's center, with the longest of them representing projected mean velocities of 258 kms"-*^ (and 164 kms"-*^ in 
the inset). The inset in the bottom left of the panel shows a zoom into the velocity field around one of Eris massive satellite companions 
(Ms = 5 X 10^ M0) relative to its own center (indicated by a plus sign in the inset). While being accreted, the satellite is also venting 
heavy elements into the surroundings. Note that, owing to the averaging process, the metallicities and velocities plotted in this figure can 
be significantly lower than the corresponding quantities for individual gas particles. 



code Sunrise (|jQnssQnll2QQ6f ), which produces spectral 
energy distributions using the age and metahicities of 
each simulated star particle, and takes into account the 
three-dimensional effect of dust reprocessing. At this 
epoch, the simulated galaxy resembles an LBG, with a 
virial mass of Mvir = 2.4 x 10^^ and a virial radius 
of i?vir = 48 kpc. The total stellar mass within i?vir 
is = 2.1 X 10^° M©. A 2D photometric decomposi- 
tion performed on the dust-reddened rest-frame i-band 
light distribution with the Galfit program (Peng et al. 
HOOE) shows the presence of an extended stellar disk 
with radial scale length = 0.6 kpc that is vigorously 
forming stars at a rate of SFR= 18 M0yr~^. The to- 
tal gas and stellar metallicities are Zg = 0.08 and 
= 0.19^0, respectively. The galaxy's gaseous disk 
(defined by all the cold, T < 3 x 10^ K gas within 
15 comoving kpc from the center) is characterized by 
12 -h log(0/H) = 8.4. This is below the mass metallic- 
i ty relation at z = f rom the Sloan Digital Sky Survey 
(jTremonti et al. I [20041 ) but in agreement with the z ^ 3 



metallicity data on LBGs from iMannucci et al. I ()2009f ). 

The complex environment of the galaxy forming re- 
gion is clearly seen in Figure [21 which shows projected 
gas density (top panel) and gas-phase metallicity (bot- 
tom panel) in a cube 500 kpc on a side. While metal- 
enriched material is seen as far as 200 (proper) kpc (i.e. 
more than 4i?vir) from the center, within the same re- 
gion we also identify 47 self-bound dwarf satellites with 
masses Ms > 10^ M© (424 with > 10^ M©), many of 
them also forming stars and polluting their surroundings. 
The total mass of heavy elements in the gas phase within 
i^vir, 2i?vir, and 3i?vir is 1.7 X 10^ M0, 2.9 x 10^ M©, 
and 3.0 x 10^ M©, respectively. A region ~ 100 kpc in 
size is enriched to metallicities above 0.03 solar. The 
extent of the metal enriched region is consistent with re- 
cent observations o f circum-galactic metals around LBGs 
(jSteidel et al.ll2010i) . 

Figure [3] shows a projected gas metallicity and ve- 
locity map in a 500 x 500 x 10 kpc^ slice. The ar- 
rows indicate the direction and magnitude of the mass- 
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Fig. 4. — Distribution of all enriched gas in the temperature- 
density plane at 2: = 3 within the simulated volume. The color 
scale indicates the mass-weighted metallicity. 

weighted peculiar velocity field relative to the center of 
Eris. Galactic winds can be clearly seen propagating 
perpendicularly to the disk (seen edge-on in this projec- 
tion) well beyond the virial radius, with average veloci- 
ties (over the slice) that exceed 250 kms~^ (we shall see 
below that individual gas particle speeds can reach 800 
kms~^). The outflows have a bipolar distribution with 
a smaller opening angle near the base, similar to the 
observed morpholo gies of galactic winds in star- forming 
galaxies (iVeilleux. Cecil, fc Bland- Hawthornll2 00 5f ). In- 
flows along large-scale fllaments bring in both pristine 
gas (gray arrows in the flgure), as well as material pre- 
enriched by dwarf satellite companions (colored arrows 
on the left side of the flgure and in the inset). 

It is instructive at this stage to look at the mass- 
weighted distribution of all (inside and outside the main 
host) enriched gas at z = 3 in the temperature-density 
plane. Figure |4] indicates that metals are spread over a 
large range of phases, from cold star- forming material at 
T < 3 X 10^ K and n > usf = 5 atoms cm~^ (corre- 
sponding to S = p/pmean > 3 X 10^ at 2: = 3) to hot 
T > 10^ K low density S = 3 intergalactic gas that can- 
not cool radiatively over a Hubble time. The black strip 
in the lower left corner of the flgure marks the pristine, 
adiabatically cooling IGM, while the colored swath in 
the lower right corner shows dense, metal-rich gas in the 
galaxy disk cooling down below 10^ K. Hot enriched gas 
vented out in the halo by the cumulative effect of SN 
explosions can be seen cooling and raining back onto the 
disk in a galactic fountain. Intergalactic gas in the range 
1 < 6 < 10 shows a strong positive gradient of metallicity 
with increasing temperature, and has the largest range 
of metallicities, extending from solar all the way down to 
zero. 

A census of all the gas-phase metals in the cold (T < 
3 X 10"^ K), warm (3 x lO'^K < T < 3 x 10^ K), and hot 
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Fig. 5. — Top panel: Evolution with redshift of the mass fraction 
of heavy elements in the gas phase within the simulated volume. 
The three curves shows the metal fraction of cold (T < 3 x 10"^ 
K), warm (3 x lO^K < T < 3 x 10^ K), and hot (T > 3 x 10^) 
interstellar and circumgalactic gas. Bottom panel: Cumulative 
metal mass fraction in each gas phase at z = 3 as a function of 
overdensity. 

(T > 3 X 10^ K) interstellar and circumgalactic medium 
within the simulated volume is depicted in Figure [5] (top 
panel) as a function of redshift. The metal fraction of 
cold gas drops from about 70% at redshift 8 to 50% at 
z = 3. Hot gas is the second most important reservoir 
of gas-phase heavy elements, with a metal fraction that 
increases from 15% at z = 8 to 40% at z = 3. This phase 
would remain undetected in UV spectroscopic studies of 
high redshift galaxies, and would therefore contribute t o 
the "missing; metals" CPettinil l2006l : iBouche et"ani2QQ7[ ) . 
The metal mass fraction in warm gas remains around 
10% at all epochs. 

More than 40% of all gas-phase metals at z = 3 lie 
outside the virial radius: while cold metal-rich material 
traces large overdensities within the main host, about 
50% of all warm and 70% of all hot metals are found 
in low density (5 < 30 regions beyond the virial radius, 
a point illustrated in the bottom panel of Figure [5l In- 
tergalactic metals are characterized by a strong temper- 
ature gradient with overdensity, as the metal-weighted 
temperature climbs from 10^ K at J = 1 to 2 x 10^ K at 
5 = 10. 

The metallicities of the individual cold, warm, and hot 
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Fig. 6. — Mean metallicity of cold, warm, and hot gas in the simulated volume as a function of gas overdensity in four redshift bins. The 
thin blue line shows the total gas metallicity. 

bution at z = 3 shows a positive density gradient above 
an overdensity of log 5 = 2, with little redshift evolu- 
tion as heavy elements removed from high density regions 
are steadily replenished. The mean total metallicity ex- 
hibits a plateau around log J = 1 — 2, where the metals 
transported by SN-driven winds accumulate, and drops 
quickly below log Z/Zq = — 2 at overdensities below a 
few. A low-density peak is more clearly seen in the metal- 
Hcity distribution of cold, warm, and hot gas: in the case 
the warm and hot medium the peak is shifted towards 
underdense 5^1 regions. At z = 3, the mean metal- 
licity of cold gas drops quickly below log Z/Zq = — 2 at 
overdensities ^ < 10. 

While our "zoom-in" simulation focuses on the circum- 
galactic medium of a single massive system and cannot 
describe the enrichment history of the IGM as a whole 
by a population of galaxies with different masses and 
star formation histories, it is still instructive to briefly 
compare our results with those of very recent, much 
lower resolution, large cosmological volume simulations 
of IGM pollution that use different su b-g:rid prescrip- 
tions for generating galac t ic outflows (IWiersnia et all 
2QQ9I: ICenfc Chisaril [2011 iQppenheimer et al. I l2Qllf r 
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Fig. 7. — Fractional amount of metals at redshift 3 at over- 
density 6(z = 3) that was added to gas particles before redshift 
z = 4,5,6.1,6.8,7.9. 



gas components, as well as that of the total gas phase, 
are shown in Figure [6] as a function of gas overdensity 
in four redshift bins. The gas-phase metallicity distri- 



Cen fc Chisaril ()2Qlll ) inject thermal energy and met- 
als from SN feedback on tens of kiloparsecs scales, do 
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Fig. 8. — Projected gas metallicity map of Eris circum-galactic medium at 2; = 3 in a 500 x 300 x 300 kpc^ (proper) box. The galaxy center 
and its virial radius are indicated by the plus sign and the black circle, respectively. The color coding indicates the mean gas metallicity in 
a column of area dxdy = 7.1 x 4.3 kpc^. The different panels highlight different enrichment redshift ranges. Gas enriched at earlier epochs 
generally lies at larger distance to the main host's center, as galactic outflows transport metals from the dense regions of star formation into 
the circum-galactic medium on cosmological timescales. At galactocentric distances > 100 kpc the IGM is enriched mostly by satellites. 
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Fig. 9. — Gaseous metal mass in 2.5 kpc thick radial shells at 
varying distances from the center at z = 3. These metals were 
released by z = 3 from Eris' main host (solid lines) and its satellites 
galaxies (dashed lines). The colors indicate metals produced at 
different enrichment epochs. 

not turn off hydrodynamic coupling between the ejected 
metals and the surrounding gas, and, like in Eris, find 
that a large fraction of all the heavy elements in the 
gas phase at redshift 3 have temperatures in excess of 
3 X 10^ K. Their metallicity-density relation for cold gas 
also shows a low-density peak, albeit shifted towards 
underdense regions compared to Eris. For comparison. 



in the momentum-conserve d wind imp lementation of ki- 
netic feedback by Oppenhei mer et al. I ([2011) (where hy- 
drodynamic forces are temporarily turned off) there is no 
low-density peak in the metallicity-density relation and 
metals largely reside in cool gas. IGM metals reside pri- 
marily in a warm-ho t component in the simulations of 
IWiersma et al.l (j2QQ9f ) , who also use kinetic feedback but 
with non-decoupled wind models. 

4. THE ORIGIN OF CIRCUM-GALACTIC METALS 

In the absence of metal diffusion (jShen et al.l l2QlQf ). 
the SPH technique allows us to trace back in time the 
enrichment history of every gas particle. In this and 
the following section we use simulation outputs from 
the ErisMC run and the Amiga's Halo Finder (AHF, 
iKnollmann fc K nebe'2QQ9^ to identify the age of the met- 
als observed at redshift 3 and the sources of pollution - 
whether dwarf satellite companions, dwarf progenitors, 
or the main host. 

Figure [3 shows the fraction of metals at redshift 3 that 
was released to gas particles at overdensity S{z = 3) be- 
fore redshifts z = 4, 5, 6.1, 6.8, 7.9. The epoch at which a 
gas element is enriched is clearly a sensitive function of its 
overdensity at some later time. The diffuse IGM is typ- 
ically enriched earlier than high density reg:ions, a trend 
that i s in agreement with the r e sults of Iwiersma et al.l 
(|2QlQr ) and lOppenheimer et all (j2Qllf ). More than 50% 
(35%) of all 2; = 3 metals at the average density were 
synthesized before z = 5 (z = 6), while newly produced 
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Fig. 10. — Metal mass in Eris' circum-galactic medium as a function of the mean source halo mass (Men)- Each panel indicates the 
physical distance of the gas from the main host's center at 2; = 3. Solid lines: main host plus satellites. Dashed lines: main host only. The 
mass bin in all panels is A (log (Men)) = 0.1. 



metals are mostly confined to high over densities. 

There are two possible causes for this "outside-in" (a 
term we borrow from Oppenheim er et al. , 2011) enrich- 
ment of Eris' circum-galactic medium: 1) metals released 
at earlier times into the high-density regions of Eris main 
host and transported into the IGM via galactic winds on 
timescales that are comparable to the age of the universe 
at z = 3; and 2) metals released at earlier times into the 
high-density regions of Eris' dwarf satellites and shed 
into the surrounding intergalactic and circum-galactic 
medium before and/or during infall. Eris' satellite galax- 
ies are referred as as "satellite companions" if they are 
still orbiting outside i?vir at z = 3, and "satellite progen- 
itors" if they have been accreted by the main host before 
redshift 3. To identify the time, location, and source 
of enrichment of a given gas particle, we de f ine a metal 
mass- weighted redshift as in IWiersma et al.l ()2Q1Q| ) : 



(^en) 



(4) 



where Amz,i is the metal mass gained by the gas particle 
in an enrichment event at redshift > 3, and Amz,i 
is the total metal mass of the particle at 2; = 3. For every 
enrichment event we also derive the mass of the (satellite) 
halo in which the gas particle resides (i.e. the particle 
belongs to the AMIGA group of the satellite halo), Mh^i 
and the distance between the gas particle and the center 
of Eris, di. Analogously to equation (j4]), we then define a 
metal mass- weighted source halo mass and distance from 



Eris center as 



and 



(Men) = 



{Den) = 



Amz,iMh,i 
Amz,i 

Amz,idi 
E^^rnz,^ ' 



(5) 



(6) 



respectively. We have separated gas particles into differ- 
ent groups according to their (zq^) and plotted in Figure 
[8] the projected metallicity of each group at redshift 3. 
Metals within Eris' virial radius are clearly "younger", 
i.e. they are characterized by an enrichment redshift 
(zen) between 3 and 3.5. Because of the long wind prop- 
agation time, "older" metals with (zen) between 4 and 
5 are spread over 100 (physical) kpc perpendicularly to 
Eris' disk. Low-metallicity gas in this enrichment red- 
shift range can be seen as far as 200 kpc from the main 
host center as it is ejected from from satellite compan- 
ions (see also Figs. [2] and [3]). There is little material 
contaminated by metals at (^^en) > 5 within Eris' virial 
radius. Late (zen < 5) galactic "superwinds" - the result 
of recent star formation in Eris' main host - are found 
to account for less than 9% of all the metals observed 
beyond 2i?vir , 

Figure [9] sheds light on the role played of satellites 
(both progenitors and companions) in contaminating 
Eris' circum-galactic medium. It shows the total mass of 
heavy elements released by the main host and its satel- 
lites as a function of distance from Eris'center at redshift 
3. About 60% of all the metals within 100 kpc of the 
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Fig. 11. — Metal mass distribution at z = 3 as a function of the comoving mean enrichment distance {-Den}- The dotted vertical lines 
and the legend in each panel mark the comoving distance interval from Eris' center of the enriched gas under consideration at z = 3. Solid 
lines: main host plus satellites. Dashed lines: main host only. 



center originate from the main host, and the rest from 
its satehites. Beyond 100 kpc, the satehites start domi- 
nating the metal budget. Both the host and the satehites 
contribute to the recent {z < 4) pohution of gas within 
the virial radius. Older metals within i^vir typically form 
in satellite progenitors, collect along the filaments into 
the main host, and are not blown away by galactic out- 
flows. Note how, within 150 kpc or so, the distribution of 
metals from satellites is rather smooth and follows that 
from the main host, an indication that gas polluted by 
star formation in satellite progenitors is stirred up and 
well mixed with Eris' galactic outflows. Spikes due to 
individual satellite companions can be seen beyond 85 
kpc. 

It is interesting at this stage to look at the masses of the 
satellites that contribute to the enrichment of the circum- 
galactic medium. Figure [10] shows the metal- weighted 
mean halo mass, (Men) (defined in eq. [5]), for gas at 
different physical distances from Eris' center. Most of 
the satellites' metals come from systems more massive 
than 10^ M0, with the peak of the distribution typically 
around 10^-^ M©. Satellites smaller than 10^-^ M© do 
not cause significant pollution as they are unable to form 
many stars because of feedback. Also, halos smaller than 
10^ M0 are resolved by less than 1,000 particles in our 
simulation, and the inability to properly resolve high den- 
sity star- forming regions in these objects may result in a 
numerical suppression of star formation (as our recipe is 
based on a high gas density threshold). The right bottom 
panel shows how gas beyond 150 kpc is enriched only by 
dwarf galaxy companions. Figure [11] provides some in- 



sight on the overall transport of enriched gas by showing 
the metal mass distribution as a function of the enrich- 
ment distance (as defined in eq. [6]). We use comoving 
distances to highlight departures from pure Hubble flow. 
Most host metals are released within 50-100 comoving 
kpc from the center, with those found beyond the virial 
radius at 2; = 3 originating earlier from strong galactic 
winds launched closer to the center. Metals from satel- 
lites can have very different kinematics. Some are ejected 
directly into the IGM by satellite companions, other, ini- 
tially deposited in the halo of the main host by infalling 
dwarf progenitors, become part of Eris' galactic outflows. 
We find that about 40-50% of all the satellite-produced 
metals found in a given distance intervals at z = 3 were 
produced at larger distances and then transported in- 
wards to their current location. For example, 42% of all 
satellite metals found at z = 3 in the range 400-600 (co- 
moving) kpc from Eris' center have a mean enrichment 
distance larger than 600 kpc. All metals found beyond 
600 comoving (150 physical) kpc were ejected by satellite 
companions that have not been yet accreted by the host. 

To better grasp the kinematics of the material enriched 
by satellite progenitors (i.e. dwarf galaxies that have 
all been accreted by the main host before z = 3), we 
plot in Figure [12] a metal column density map of Eris' 
circum-galactic environment at three different observer 
redshifts, Zobs = 3,4,5. All the heavy elements shown in 
the figure were produced by satellite progenitors at epochs 
5 < z < 7. All distances are comoving and the metals 
are separated according to their radial peculiar velocities 
relative to the center of the main host, Vr^ in inflowing 
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Fig. 12. — Metal column density map of Eris' circum-galactic environment at three different observer redshifts, 2;obs = 3,4,5. All the 
heavy elements shown in the figure are produced by satellite progenitors at epochs 5 < z < 7. The projected region has a comoving size of 
2 X 1.2 X 1.2 Mpc^ and is centered on the main host. The center of the host galaxy and its virial radius are indicated by the plus sign 
and black circle, respectively. The color scale indicates the logarithmic of total metal mass a column dxdy = 28.6 x 17.1 kpc^ comoving, 
in units of Mq. Left panels: inflowing material. Right panels: outflowing material. 



{vr < 0, left panels) and outflowing {vr > 0, right pan- 
els). The projected mass distribution is shown in the 
gray scale. Inflowing metals at Zobs = 5 have contam- 
inated the vicinity of their satellite hosts, and are now 
being accreted onto the main host along the fllamentary 
structure. They are subsequently dispersed during the 
infall and tidal disruption of their satellite hosts, become 
entrained in the galactic wind of the main host, and are 
ultimately ejected in a bipolar outflow. 

5. PROPERTIES OF GALACTIC OUTFLOWS 

In this section, we study the properties of Eris' galactic 
onflows and compare them with the observations as well 
as with other galactic wind models adopted in cosmolog- 
ical simulations. 

5.1. Outflow velocity 



An example of the evolution with time of the ra- 
dial velocity of polluted material is shown in Figure 
fT3l where we have selected gas particles that were 
enriched for the last time in the main host at red- 
shift 5, and that are unbound at 2; = 3. Gas ac- 
cretes with negative peculiar velocity onto the host, 
is enriched and accelerated to outflow velocities of 
a few hundred kms~^ by blastwave feedback, and 
slows down as it sweeps the ISM and halo material. 
The peak velocity is comparable to those observed 
in high redshift LBGs (e,g., Adelberger et al. 2003* 
Veilleux. Cecil, fc Bland- Hawtho rn 2005; Steidel e t al.1 
2010^ It is also comparable to the veloci- 
ties adopted in various kinetic feedback models 
(e.g., 'Sprinffel & Hernauist' '2003'; 'Qppenhei mer fc Pavel 
[2006: Daha Vecchia & Schaye 2008; Wiersma ~al.ll2009l : 
iChoi fc Nagaminel 12OII1) . After reaching a peak value. 
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Fig. 13. — Evolution of the radial velocity of gas particles last 
enriched at redshift 5 in the main host. Peculiar and Hubble flow 
velocities are indicated with the black and red curves, respectively. 
The solid, dotted, and dashed lines show the median, 5 percentile 
and 95 percentile values, respectively. 



Fig. 14. — Peak radial velocity of gas particles that become un- 
bound after they are enriched in the main host, as a function 
of their enrichment redshift. Each dot represents a gas particle. 
The black and blue solid lines show the mass-weighted and metal- 
weighted average values, respectively. The colored dashed lines 
show the standard deviation from the mean. 



the gas' proper velocity declines rapidly, and particles 
are swept by the Hubble flow as they move farther from 
the center into the IGM. This behaviour is typical of all 
enriched gas particles. The rapid decline of the proper 
velocity is also seen in the kinetic feedback model of 
iDalla Vecchia fc S chave (2008), where outflowing parti- 
cles are allowed to interact hydrodyamically with the 
ISM, and pressure forces within the disk significantly de- 
crease the wind speed. 

Observations of local starburst g alaxies (e.g., 
[Schwartz & Martin 2004; "Rupke et alj I2005D have 
shown that outflow speeds are correlated with halo 
masses and star formation rates, and that the correla- 
tion flattens out for SFR> 10 M© yr"^ (|Rupke et al.l 
l2005l ). To investigate the existence of such relationship 
in the Eris simulation, we plot in the upper panel of 
Figure [M] the peak velocity of gas particles that become 
unbound after they are enriched in the main host, as 
a function of their enrichment redshift. We choose the 
gas peak velocity as this is close to the flow speed when 
the wind is launched, and is relatively unaffected by 
interactions between the outflow and the ISM/gaseous 
halo. Note that, because of the limited time resolution, 
we may actually underestimate the peak velocity of 
gas particles that reach a maximum speed and then 
slow down significantly within the time interval of 
two simulation outputs {At > 140 Myr). The scatter 
plot shows that enriched gas particles can often reach 
velocities in excess of 600 kms~^, with a few of them 
moving as fast as ~ 800 — 1000 kms~^, a value that is 
consistent with the higs:h est velocity n i ateria l observed in 
LBGs at 2 < z < 3 by ISteidel et all (j2010f ). The mass 
and metal-weighted outflow average speed typically 
ranges between 200 and 400 kms~^. From redshift 9.3 
to 3.0, the total mass of the main host halo increases 
from 5.0 x 10^ M© to 2.4 x 10^^ M©. The mean outflow 
speed increases for 5 ^ z < 9.3 and decreases again 
at z < 5, so there is no obvious correlation between 
halo mass and outflow speed in Eris. Similarly, we find 
only a weak correlation between wind velocity and star 
formation rate. 
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Fig. 15. — Mean metallicity of inflowing (dashed lines) and out- 
flowing (solid lines) material as a function of redshift. The average 
is taken over all gas particles within a thin shell of radius R^ir cind 
thickness 0.02i?vir- 

5.2. Metallicity 

Figure [15] shows the metalHcity of inflowing and out- 
flowing material at Eris' virial radius as a function of red- 
shift. Galactic outflows are enriched to a typical metal- 
licity in the range 0.1 — 0.2 Zq since redshift ~ 9, with 
little dependence on cosmic time. The mean metallicity 
of inflowing gas shows a more marked evolution, as it 
increases by about one dex in the interval 9 ^ z ^ 3.5. 
This is a consequence of more processed material falling 
back onto Eris in "galactic fountains" as well as being 
accreted via infalling satellites. The gas metallicity of 
inflowing gas atz<5(Z~0.01 Zq) is typical of the 
metallicity ob served in Damped Lya and Lyman-Limit 
systems f e.g.. I Wolfe et al.ll20Q5[ ). 

5.3. Mass-loading factor 

The mass loading factor characterizes the amount 
of material involved in a galactic outflow, and is 
defined as = M^/SFR, where is the rate 

at which mass is ejected. Observations of galac- 
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Fig. 16. — Mass loading factor as a function of cosmic time or 
redshift. The mass flux was measured at the virial radius (black 
solid and dashed curves), half the virial radius (blue curves) and 
one flfths of the virial radius (red curves). The solid and dashed 
lines indicate the total mass loading factor versus the loading factor 
of enriched gas only. 

tic outflows powered by starbursts suggest a wide 
range of mass loading factors, 77 = 0.01 — 10, 
with no obvious co rrelation with the star formation 
rates of their hosts (jVeilleux. Cecil, fc Bland-HawthornI 
[2005i) . In low-resolution cosmological simulations, the 
mass loading factor is one of the input parameters 
(e.g., Springel & Hernauist 2003; Oppenheimer & Dave 
[2OO6 : Dalla Vecchia & Schaye 2008; Choi & Nagamine 
l201l[ ). Most models adopt a constant mass load- 
ing factor, except for the momentumr-driven model 
(jOppenheimer fc Pavel l2006l . l2009l ). where it is assumed 
to be inversely proportional to the host galaxy veloc- 
ity dispersion a, 77 = ao/cr. In our SN-driven blastwave 
feedback scheme there is no specific parameter for mass 
loading. In this section we explore the mass loading fac- 
tor in Eris and its variation with star formation and halo 
mass. 

We compute the wind mass loading at a distance r 
from the center as 

1 ^ 

Mw{r) = —y^miVr^i (7) 

i 

where Ar is the thickness of the radial shell, N is the 
total number of outflow gas particles in the shell, rrii is 
the mass of particle i and Vr^i its outflow radial velocity 
(relative to the host's center). We have measured the 
mass flux at half the virial radius and at the virial ra- 
dius, and divided it by the instantaneous SFR ignoring 
any time delay between star formation and large-scale 
outflow. Figure [16] shows that the mass loading factor 
stays between 0.1 and 1.2 at all redshifts. These val- 
ues are smaller than those assumed in kinetic feedback 
and in momentum-driven wind models. We find no clear 
correlation between mass loading factor and the mass or 
velocity dispersion of the host galaxy. 

6. SUMMARY 

We have presented a detailed study of the metal- 
enriched circum-galactic medium of a massive galaxy at 
z = 3 using results from the Eris suite of new cosmo- 
logical hydrodynamic "zoom-in" simulations in which a 



close analog of a Milky Way disk galaxy arises at the 
present epoch. Our approach to understanding the role 
of inflows, star formation feedback, and outflows in gov- 
erning the gaseous and metal content of galaxies and 
their environment, is different to that of many recent 
theoretical efforts: Eris' extreme mass and spatial reso- 
lution allows us to follow self-consistently the venting of 
metals by small progenitor dwarf satellites and the trans- 
port of heavy elements from their production sites into 
the environment. The reference run adopts a blastwave 
scheme for supernova feedback that generates galactic 
outflows without explicit wind particles, and a star for- 
mation recipe based on a high gas density threshold. 

We have found that Eris' metal-enriched circum- 
galactic medium extends as far as 4 virial radii (about 
200 physical kpc) from its center. Approximately 41%, 
9%, and 50% of all gas-phase heavy elements are hot 
(T > 3 X 10^ K), warm (3 x lO^K > T > 3 x 10"^ K), 
and cold (T < 3 x 10"^ K), respectively. More than 40% 
of all gas-phase metals lie outside the virial radius: while 
cold metal-rich material traces large overdensities within 
the main host, about 50% of all warm and 70% of all hot 
metals are found in low density ^ < 30 regions beyond 
Eris' virial radius. Intergalactic metals are characterized 
by a strong temperature gradient with overdensity, as 
the metal- weighted temperature climbs from 10^ K at 
S = 1 to above 2xlO^Kat^ = 10. SN-driven winds are 
able to transport metals to regions at J < 10 where they 
accumulate, creating a peak in the cold gas metallicity- 
density relation at Z/Zq ^ 0.05. The mean metallicity 
of cold gas drops quickly below log Z/Zq = — 2 at over- 
densities 6 < 10. 

We have identified three sources of heavy elements in 
the simulated region: 1) the main host, responsible for 
60% percent of all the metals found within 3i?vir , and for 
none of those found beyond 3i?vir ; 2) its satellite progen- 
itors - systems accreted by the main host before redshift 
3, which shed their metals before and during infall and 
are responsible for 28% of all the metals within 3i?vir ; 
and its satellite dwarf companions - still orbiting outside 
the virial radius at redshift 3, which give origin to 12% 
of all the metals within 3i?vir and 95% of those beyond 
3i?vir • Late {z < 5) galactic "superwinds" - the result of 
recent star formation in Eris - account for only 9% of all 
the metals observed beyond 2i^vir , the rest having been 
released at redshifts 5 ^ z < 8. Th ese findings confirm 
the ideas put forward bv lPorciani fc Mad au (2005), that 
"pregalactic" enrichment from dwarf satellite systems 
may contribute significantly to (and dominate at large 
distances) the pollution of the circum-galactic medium 
around LBGs. 

Substantial amounts of heavy elements are generated 
at a larger distance from the main host's center than 
their current location, and subsequently are accreted 
by the host along filaments via low-metallicity cold in- 
flows. Galactic outflows have velocities of a few hun- 
dred kms~^. The outflows decelerate rapidly, and the 
resulting long metal transport timescales produce an age 
gradient in metals as a function of distance from the 
main host. The outflow mass- loading factor as a func- 
tion of redshift ranges between 0.1 and 1.2, and shows 
no correlation with the mass of the host over the interval 
5.0 X 10^ M« < Mvir < 2.4 X 10^^ Mc:^. 
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As a Lagrangian particle method, SPH does not in- 
clude any implicit diffusion of scalar quantities such 
as metals. In the absence of some implementation 
of diffusion, metals are locked into specific particles 
and their distribution may be artificially inhomogeneous 
(|Aguirr e et al. 2005; Wiersma et al. 2009). A number of 
simulations of even higher resolution than Eris and in- 
cluding a scheme for turbulent mixing that redistributes 
heavy elements and thermal energy between the outflow- 



ing material and the ambient gas (a la lShen et al.ll20fQh 
are in the making. 
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